APPENDIX A 



1 package ihn.dialex; 

2 import ihn.merkato.*; 

3 import ihn.math.Gaussian; 

4 import ihn. math. Normal; 

5 import java.uti I. Hashtable; 

6 import java.util. Random; 

7 /** 

8 *@version SRevision: 1.6 $, $Date: 1998/12/22 05:21:22 $ 

9 *@author Nemo Semret 

10 ♦/ 

n public class ReservationCalculator extends Hashtable { 

12 intV2=l; 

13 intC,B; 

14 float lambda. mu,r; 

15 intdt; 

16 int horizon; 

17 public ReservationCalculator(int C.int B,float mu, float lambda, 

18 float r, int horizon, int dt) { 

19 this.C = C; 

20 this.B=B; 

2 1 this.mu-mu; 

22 this.lambda=lambda; 

23 this.r=r; 

24 this.horizon= horizon; 

25 this.dt = dt; 

26 System.out.println(this); 

27 » 
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public float interArrivalO { 
return (float) (1.0/lambda); 

I 

public String toString() { 
return super.toString()+" C = "+C +" B="+B +" mu="+mu + 

" lambda="+lambda+" r="+r+" horiz="+horizon+" dt="+dt; 

> 

float A,K3,D,K4; 
/•* 

♦Compute beta, K.1, K2 for each m 

*and A,K3,D,K4. 

*/ 

public void init() throws ArithmeticException { 

float [] Kl; 
float [] K2; 
float [] beta; 



Kl =new float [B+l]; 

K2 = new float [B+l]; 

beta = new float [B+l]; 

float [] wO = new float [B+l]; 

float [] wl = new float [B+l]; 

Gaussian phi = new Gaussian(0,l); 

float root2overpi = (float)Math.sqrt(2.0/Math.PI); 

float rho = lambda/mu; 

float zob = rho/(l+rho); 

float plop = l/(l+rho); 

A- K3- D= K4=0; 
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55 float rC = r*C; 

56 float rootrC = (float)Math.sqrt(rC); 

57 float u; 

58 

59 for (int m=0 ; m<=B; m++, p!op=plop*zob) { 

60 if(m>=4*V2) { 

61 u = (1- (float)Math.sqrt(l-(float)(4*V2)/m))/2; 

62 wO[m]=(float)Math.sqrt(u*m/(l-u)); 

63 wl[m]=(float)Math.sqrt((l-u)*m/u); 

64 Kl [m]= 2*phi.integral(w0[m],wl tm],(float)l .0); 

65 K2[m]= root2overpi*((float)Mathxxp(-wO[m]*wO[rn]/2)*(wO[rn]/m - (float) 1 .0/w0[m]) - 

66 (float)Math.exp(-wl[m]*wl[m]/2)*(wl[m]/m - (float) 1.0/wl[m]))+2*Kl[m]; 

67 if(m<B) beta[m]=plop; 

68 else beta[m]= plop*rho; 

69 A -=beta[m]*rC/m; 

70 K3 +=beta[m]*rC*((m-l )*K1 [m]/m + K2[m]/m); 

7 1 D+= beta[m]*rootrC/m; 

72 K4 -=beta [m]*Kl [m]*rootrC; 

73 } 

74 } 

75 System.out.print]n("A="+A+" K3="+K3+" D= M +D+" K4="+K4); 

76 » 

77 /** 

78 * Generates a forward sample path of the process P. 

79 * @param pO The initial price P[0]=po; 

80 * @param T The duration. 
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81 * @retum An array of size l+floor(T/dt) 

82 */ 

83 float [] fwdSamplePath(float pO, int T, int dt ) { 

84 float dP=0, dW=0; 

85 float []P = newfloat[l+T/dt]; 

86 // Normal norm = new Normal(l 00); 

87 Random rand = new Random(System.currentTimeMillisO); 

88 P[0]=p0; 

89 int i,t; 

90 for( t=dt, i=l ; t<=T; t+=dt) { 

91 // dP = A*P[i-l]+K3*dt + (D*P[i-l]+K4)*norm.nextSample()*(float)Math.sqrt(dt); 

92 dP = A*P[i-l ]+K3*dt + (D*P[i-l ]+K4) 

93 *(float) (rand.nextDouble() <0.5 ? -1 .0 : 1 .0) 

94 *(float)Math.sqrt(dt); 

95 P[i]=P[i-l]+dP; 

96 } 

97 return P; 

98 > 

99 /** 

100 * For each time t between the current and T, we get an array 

101 * which contains the histogram of P[t]. 

102 *@param NBINS Number of bins in the histogram, 

103 *@param pO The current price. 

104 

105 •/ 

1 06 int [][] monteCarlo(float pO, int T) { 
107 

1 08 int [][] histog = new int[l+T/dt][NBINS]; 

109 float [] samplepath ; 
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110 intn.t.i; 

111 int p; 

112 for ( n=0; n<n samples; n++) { 

113 samplepath = fwdSamplePath(pO,T,dt); 

114 for (1=0, i=0; t<-T; i++, t+=dt) { 

115 P=(int)( samplepath[i]*NBINS); 

116 p = p<0?0:p; 

117 p= p>NBINS-l? NBINS-1 :p; 

118 histog[i][p]++; 

119 > 

120 } 

121 return histog; 

122 ) 

123 /» 

124 * Number of bins in each histogram 

125 */ 

126 intNBINS=10; 

127 /« 

128 * Number of sample paths for each fastFwd 

129 •/ 

130 int nsamples=100; 

131 public synchronized void fastFwd(float pO) { 

132 if( get(new Integer((int)(pO*NBINS))) = null) { 

133 System.outprint("FFing "+hori2on+" for price "+p0); 

134 long a=System.currentTimeMillis(); 

135 put(new Integer((int)(pO*NBINS)), monteCarlo(pO, horizon)); 

136 System.out.println( ,, done in "+(System.currentTimeMillis()-a)); 

137 } 
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138 



139 /** 

1 40 *@param pO The current market price. 

141 *@param b The reservation for which the fee is computed. 

142 •/ 

143 public Reservation fee(Bid b, float pO) { 

144 intx,t,i; 

145 float fee=0; 

146 float duration = Math.min(b.qty,horizon); 

1 47 int [][] histogram = (int [][]) get(new Integer((int)(pO*NBINS))); 

148 if(histogram — null) { 

149 fastFwd(pO); 

150 histogram = (int [][])get(new Integer((int)(pO*NBINS))); 

151 } 

1 52 for(i=0, t=0; t<=duration; t+=dt, 

153 for( x=0; x<NBINS; x++) 

1 54 fee += (Math.max((x+0.5)/((float)NBINS) - b.price.O) 

155 *histogram[i][x])/nsamples; 

156 fee=fee*dt; 

157 Reservation foo=new Reservation(b); 

158 b.qty =duration; 

159 foo.fee= fee; 

160 // System.out.println(pO+" "+foo.price+" "+foo.qty+" "+foo.fee+" "+foo.bidderid); 

161 return foo; 

162 } 



163 > 
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